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Anomalously hindered E2 strength B{E2]l\ — > + ) in 16 C 
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The electric quadrupole transition from the first 2 + state to the ground + state in 16 C is studied 
through measurement of the lifetime by a recoil shadow method applied to inelastically scattered 
radioactive 16 C nuclei. The measured lifetime is 75 ± 23 ps, corresponding to a B(E2;2^ — > + ) 
(f) ■ value of 0.63 ± 0.19 e 2 fm 4 , or 0.26 ± 0.08 Weisskopf units. The transition strength is found to be 

anomalously small compared to the empirically predicted value. 
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PACS numbers: 23.20.Js, 21.10.Tg, 29.30.Kv 



Quadrupole strengths are fundamental quantities in 
probing the collective character of nuclei. The en- 
hancement of the electric quadrupole (E2) transition 
strength with respect to that of single proton excitation 
may reflect large fluctuation or deformation of nuclear 
charge [lj. One of the important E2 transitions in an 
even-even nucleus is that from the first 2 + (2*) state to 
the ground state (0+ s ), the reduced transition probabil- 
ity B(E2) of which has long been a basic observable in the 
extraction of the magnitude of nuclear deformation or in 
probing anomalies in the nuclear structure. With recent 
advances in techniques for supplying intense beams of un- 
stable nuclei, several exotic properties such as magicity 
loss HJ. 0- Q have been discovered in neutron-rich nuclei 
through measurements of E2 strengths. 

The present Letter reports lifetime measurements of 
the 2j state of the neutron-rich nucleus 16 C. The lifetime 
is inversely proportional to B(E2). A simple model of a 
nucleus as a quantum liquid-drop well describes the sys- 
tematic tendency that B(E2) varies in inverse proportion 
to the excitation energy E(2 J [) of the 2f state Q- For 
carbon isotopes, when N changes from the magic num- 
ber 8 to 10, E(2f) decreases dramatically from 7.012 
to 1.766 MeV : 6:]. It can then be anticipated that 16 C 
(JV = 10) would have a much larger B(E2) than that of 
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14 C (N = 8). Unexpectedly, a remarkably small B(E2) 
was found for 16 C in the present work. The observed 
value, in Weisskopf units, turns out to be far smaller 
than any other B(E2) measured on the nuclear chart. 

In the present experiment, a new technique was em- 
ployed to measure the lifetime of an excited state pop- 
ulated in inverse- kinematics reactions. The technique 
essentially followed the concept of the recoil shadow 
method (RSM) 0, in which the emission point of the 
de-excitation 7-ray is located and the 7-ray intensity is 
recorded as a function of the flight distance of the de- 
exciting nucleus. As the flight velocity of the de-exciting 
nucleus is close to half the velocity of light, the flight dis- 
tance over 100 ps corresponds to a macroscopic length of 
about 1.7 cm. Thus, the present shadow method provides 
a wide range of applicability, extending to lifetimes of as 
short as a few tens of ps. In particular, the method is 
useful for determining the B(E2) value of Z < 10 nuclei, 
at which the use of intermediate-energy Coulomb exci- 
tation 0, Q may suffer from contamination of nuclear 
excitation. 

The experiment was performed at the RIKEN accelera- 
tor research facility. A secondary 16 C beam was produced 
through projectile fragmentation of a 100-MeV/nucleon 
18 primary beam, separated by the RIPS beam line 0. 
The 16 C beam was directed at a thick, 370-mg/cm 2 9 Be 
target placed at the exit of the RIPS line for inelastically 
excitation. Particle identification for the secondary beam 
was performed event-by-event by means of the time-of- 
flight (TOF)-A-E method using a 1.0 mm-thick plastic 
scintillation counter (PL) located 180 cm upstream of 
the target. Two sets of parallel plate avalanche coun- 
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FIG. 1: Schematic of experimental setup. A beryllium target 
is surrounded by a 5 cm-thick lead wall, and two cylindri- 
cally layers of Nal(Tl) scintillators (Rl and R2) are placed 
upstream of the target. 



ters (PPACs) were also placed upstream of the target to 
record the position and angle of the projectile incident 
upon the target. The 16 C beam had a typical intensity of 
2x 10 5 particles per second and purity of about 97%. The 
average energy across the target was 34.6 MeV/nucleon. 

Particles scattered from the target were identified by 
the AE-E-TOF method using a plastic scintillator ho- 
doscope located 86 cm downstream of the target. The 
hodoscope, with an active area of 24x24 cm 2 , consisted 
of a 2.0 mm-thick AE plane and a 5.0 mm thick E plane. 
The scattering angles were measured by another PPAC 
placed 25 cm downstream of the target. 

The experimental setup for 7-ray detection is shown 
schematically in Fig. ^ In order to implement the RSM 
concept, a thick 7-ray shield was placed around the tar- 
get. The shield was a 5 cm-thick lead slab with an outer 
frame of 60x60 cm 2 and an inner hole of 6.2 cm <j>. The 
inner hole surrounded the beam tube housing the 9 Be 
target. For the sake of later discussion, the z-axis is de- 
fined as the beam direction, that is close to the flight 
direction of the de-exciting nucleus. The origin of the 
z-axis, z = 0.0 cm, was taken at the upstream surface of 
the lead slab. 

The 7-rays from the excited 16 C in-flight were detected 
by two rings of Nal(Tl) detectors labeled Rl and R2, 
composed respectively of 14 and 18 rectangular Nal(Tl) 
crystals with volume of 6 x 6 x 12 cm 3 . These rings circled 
the beam tube, with crystal centers set at polar angles of 
121° and 102°, respectively. The average distance from 
the target to Rl was 29.2 cm, and 28.7 cm for R2. 

The acceptance of the detectors was determined by the 
geometry of the Rl and R2 detectors with respect to the 
lead slab, and the emission point z. Specifically, the 7- 
rays emitted between z = 0.0 and 2.2 cm can reach the 
centers of the Nal crystals of Rl without passing through 
the lead slab (z = 0.0 and 0.7 cm for R2). Hence, the rel- 
ative efficiency of R2 with respect to Rl decreased as the 
emission point was moved further downstream. This dif- 
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FIG. 2: Energy spectrum of 7-rays measured at Rl for tar- 
get position z = 0.0 cm. Dashed line represents the transition 
of 1766 keV from the excited states of 16 C, and dash-dotted 
line denotes that for 2376 keV, as simulated using GEANT 
code [jj. The solid curve represents the best fit of the sim- 
ulated line and an exponential background. The dotted line 
denotes the quoted background. 



ference of geometrical efficiency results in a dependence 
of the 7-ray yield ratio (R1/R2) on the lifetime of the 
decaying state. 

The measurement of R1/R2 was conducted for two 
target positions, z — 0.0 and 1.0 cm. Figure shows 
a Doppler-uncorrected 7-ray energy spectrum measured 
by Rl for the target placed at z = 0.0 cm. The spectrum 
is dominated by a peak at around 1.5 MeV, representing 
the Doppler-shifted 1766 keV line of the 2+ — > 0+ s tran- 
sition in 16 C. A minor peak at around 2.0 MeV is also 
present, which may correspond to the known 4 + —> 2f 
transition in 16 C. As the intensity of the latter transition 
is only 2% of the former, the contribution of the minor 
peak is ignored in subsequent analysis. 

In deducing the yield ratio, a Gaussian shape for the 
peak and cither an exponential or polynomial shape for 
the background component were fitted to the 7-ray en- 
ergy spectrum. The R1/R2 ratios thus obtained for the 
1766-keV 7-rays were 1.06 ± 0.03 and 1.70 ± 0.06 for 
target positions of z = 0.0 and 1.0 cm, respectively. The 
errors are statistical errors only. The influence of the 
choice of functional for the background is negligible. 

In order to relate the R1/R2 ratio to lifetime, Monte 
Carlo simulation was performed to evaluate the detection 
efficiencies of Rl and R2 as a function of the position z 
of 7 emission. By folding the efficiency over the flight 
distance, the R1/R2 ratio can be obtained as a func- 
tion of mean-life r. The simulation was performed using 
GEANT code [|| , and involved the shape of the detector 
and the geometry of the experimental setup. The va- 
lidity of the simulation was tested by first performing a 
separate measurement in which a 22 Na standard source 
emitting 1275-keV 7 rays was placed at twenty-one po- 
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FIG. 3: Two solid lines represent r vs. R1/R2 curves 
obtained by Monte Carlo simulation for target positions of 
z — 0.0 and 1.0 cm. The hatched zones represent the exper- 
imentally deduced R1/R2 ratios for the two target positions. 
The dashed line denotes the adopted lifetime, and the two 
dotted lines represent the range of error. 



sitions from z = 0.0 to 2.0 cm. The R1/R2 ratios thus 
measured for the different geometries were all reproduced 
by the simulation within ± 3%. 

The simulation was then applied to the case of 7-rays 
emitted from de-exciting fragments in flight. This simula- 
tion employed experimentally obtained parameters such 
as the energy and emittance of the projectile, the angular 
spread due to inelastic scattering and multiple scatter- 
ing, and energy loss of the incoming and outgoing par- 
ticles. The angular distribution of the 7-ray emission 
was also incorporated. The distribution was estimated 
using ECIS79 code 0], employing the optical potential 
parameters determined for the 12 C + 12 C reaction [Tlj . 
The r vs. RI/R2 curves thus obtained for the two target 
positions are shown in Fig. [3] as solid lines. 

This figure compares the measured Rl /R2 ratios with 
the simulated curves, allowing the value of r to be de- 
termined. The vertical widths of the two hatched zones 
represent the ranges of R1/R2 values obtained experi- 
mentally. From the overlap between the simulated curves 
and the hatched zones, the values of r were determined 
as 92 ± 22 ps and 63 ± 17 ps, respectively, for data ob- 
tained at target positions of z = 0.0 and 1.0 cm. The 
results obtained for the two different target positions are 
in reasonable agreement. 

It is probable that the choice of optical potential may 
affect the calculated 7-ray angular distribution. This 
possibility was examined by repeating the distorted wave- 
born approximation calculation using a different set of 
parameters obtained for the 16 O + 12 C reaction [l^ . 
With the resultant simulation curves, the values of r were 
obtained as 93 ± 21 ps and 66 ± 16 ps for the two target 
positions, z = 0.0 and 1.0 cm. These values are close to 



those obtained using the original parameters of optical 
potential. The weighted average of these four values is 
taken as the final value. 

The uncertainty of t considered above, 19%, is at- 
tributable entirely to statistical error. Systematic errors 
may have arisen primarily from uncertainty of the target 
position, the accuracy of which was determined to 20%. 
Another source of systematic error, the 5% difference in 
lifetimes determined by the two different sets of optical 
potential parameters, should also be taken into account. 
The total systematic error thus estimated is about 25%. 
By considering all these uncertainties, the lifetime was 
determined to be r — 75 ± 23 ps, which corresponds to 
the B(E2) value of 0.63 ± 0.19e 2 fm 4 , or 0.26 ± 0.08 in 
Weisskopf units (W.u.). 

In Fig. 01a), the B(E2) value obtained for 16 C is com- 
pared with all the other B{E2) values known for the nu- 
clei with A < 50 0]. Nuclei with open shells tend to 
have B(E2) values of greater than 10 W.u., whereas nu- 
clei with shell closure of neutrons or protons or both, 
tend to have distinctly smaller B(E2) values. Typical 
examples of the latter category are doubly magic nuclei, 
16 and 48 Ca, for which the B{E2) values are known 
to be 3.17 and 1.58 W.u., respectively. The value of 
B{E2) = 0.26 W.u. obtained for 16 C is even smaller 
than these extreme cases by as much as an order of mag- 
nitude. 

The anomalously strong hindrance of the 16 C transi- 
tion can be also illustrated through comparison with an 
empirical formula based on a liquid-drop model pj. Ac- 
cording to the formula, the values of E(2^) and B(E2) 
are related by 

B(E2) sys = 6.47 x Z 2 A-° m E{2+)- 1 . 

The experimental B(E2) values relative to B(E2) sys , 
are plotted in Fig. Elb). As noted in Ref. [3, the 
B(E2) I B(E2) sys . ratios for open-shell nuclei mostly fall 
around 1.0, being confined between 0.5 and 2.0. Even for 
the closed-shell nuclei, the ratio remains larger than 0.20. 
Thus, the ratio of 0.03 for 16 C is exceptionally small, far 
smaller than for any other nuclei, including closed-shell 
nuclei. 

Recently a hindered B(E2) has been reported for 
136 Te ^4|. This nucleus represents a case similar to 16 C 
in the viewpoint of neighboring a doubly magic nucleus. 
The B(E2) value observed for 136 Te is 206 e 2 fm 4 , which 
corresponds to B(E2)/ B(E2) sys . = 0.21. Thus, the E2 
strength is indeed hindered. However, the degree of the 
hindrance is comparable to that of singly closed nuclei, 
and is far more moderate than for the case of 16 C. 

In considering the nature of the 16 C transition, it is 
instructive to invoke a simple seniority-two empirical 
model [l4| to describe the wave function of the 2\ state 
of 16 C. The nucleus 16 C has two missing protons and two 
extra neutrons with respect to 16 O. Then, the mixing of 
7r~ 2 and v 2 configurations with a residual interaction of 
V can be considered by taking the observed 2^ states of 
14 C and 18 as the basis states. The two singly closed 
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FIG. 4: (a) B(E2) values for A < 50 W.u. (b) Ratios be- 
tween B(E2) values and values calculated using the empiri- 
cal formula [3- Closed circles denote the values of 16 C, and 
open circles represent data for other open-shell nuclei. Open 
squares denote proton closed-shell nuclei, and diamonds in- 
dicate neutron closed-shell nuclei. Full diamonds represent 
double magic nuclei. 



isotopes, 14 C with N = 8 and 18 with Z = 8, have 
two missing protons and two extra neutrons with respect 
to 16 0, respectively. Thus, the 2^ state of 14 C may be 
regarded as primarily due to proton excitations, while 
that of 18 O can be attributed almost entirely to neutron 
excitations. The E(2f) values are 7.01 and 1.98 MeV, 
and the B(E2) values are 3.70 and 9.53 e 2 fm 4 (or 1.84 
and 3.40 W.u.) for 14 C and ls O, respectively. The very 
large E(2± ) for 14 C indicates a large energy gap between 
the two proton orbitals, 1^3/2 and lpi/2, while the mod- 
erate energy for 18 O may reflect unoccupied sd-shell or- 
bitals available for neutron excitation. Correspondingly, 



B(E2) for 14 C is close to 1 W.u., whereas that for 18 is 
somewhat larger. 

By applying the empirical model, \V\ = 1067 keV is 
required to generate the 1766-keV state in 16 C. The very 
high energy of 1 7r 2 } results in a configuration dominated 
byk 2 ): 

|2+; 16 C) =0.20|tt- 2 } -0.98|^ 2 }. 

Here, the interaction V between proton-holes and 
neutron-particles is repulsive in the present case |15| . 
leading to destructive interference between the two E2 
components. The calculated B{E2) for 16 C is 7.0 e 2 fm 4 , 
which falls between these components of 14 C and 18 O, 
and overestimates the experimental B(E2) value of 16 C 
by almost an order of magnitude. This is in contrast 
to the case for 136 Te, which when treated by seniority- 
two scheme provides a fairly reasonable estimate of the 
experimental B(E2) |14| . The unprecedented hindrance 
for 16 C is suggestive of an anomalous nuclear structure. 

In this regard, several unique features of the nuclei 
in the vicinity of 16 C can be noted. For example, it 
has been suggested [l6| that the energy gaps between 
single-particle orbitals such as 1^3/2 and lpi/2, may vary 
appreciably as the isotope changes from 14 C to 16 C. It 
has also been indicated based on electric quadrupole mo- 
ment measurement of neighboring nuclei 15 > 17 B 0, 0] 
that the E2 effective charges can be markedly reduced for 
extremely neutron-rich nuclei. Furthermore, delicate in- 
terplay between clustering and deformation may prevail 
in these light nuclei [l^. Elaborate calculations involv- 
ing such exotic features can be expected to be helpful in 
accounting for the hindrance of B(E2) observed for 16 C. 
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